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Ecosystems are commonly conceptualized as networks of interacting species. However, parti¬ 
tioning natural diversity of organisms into discrete units is notoriously problematic, and mounting 
experimental evidence raises the intriguing question whether this perspective is appropriate for the 
microbial world. Here, an alternative formalism is proposed that does not require postulating the 
existence of species as fundamental ecological variables, and provides a naturally hierarchical de¬ 
scription of community dynamics. This formalism allows approaching the “species problem” from 
the opposite direction. While the classical models treat a world of imperfectly clustered organism 
types as a perturbation around well-clustered “species”, the presented approach allows gradually 
adding structure to a fully disordered background. The relevance of this theoretical construct for 
describing highly diverse natural ecosystems is discussed. 
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Although the basic unit participating in ecologi¬ 
cal interactions is an individual organism, constructing 
tractable theoretical models usually requires clustering 
individuals into discrete groups within which organisms 
are treated as identical; such classification can be based, 
for example, on taxonomy, development stage, or phe¬ 
notype [1]. The resulting tension is a long-standing is¬ 
sue in community ecology: ever since Darwin [2], the 
difficulty of drawing sharp boundaries partitioning nat¬ 
ural diversity into discrete categories [3] and the real¬ 
ization that the individual-level variation can be an im¬ 
portant actor in ecological phenomena gnu made the 
“partitioned community” assumption of well-delimited, 
uniform groups highly problematic [6Hn]. 

The urgency of this issue has been highlighted over the 
last decade by the studies of microbial diversity in natu¬ 
ral environments [IUTT]. For microbes, the partitioning 
problem is intensified by the prevalence of asexual repro¬ 
duction and horizontal gene transfer [181 IE] ? i^he fact 
that even genetically identical cells can assume different 
ecological roles, e.g. in a biofilm [20]. To cope with these 
issues, considerable attention has been devoted to iden¬ 
tifying questions that can be asked and answered with¬ 
out explicitly specifying a partitioning m, e.g. solely 
in terms of community metagenome [IIl[2ll[22]. These 
approaches adopt a convenient coarse-grained viewpoint 
that is appropriate for certain functional questions, e.g. 
comparing communities [21]. However, their ability 
to describe community dynamics is necessarily limited. 
For that purpose, the partitioned community assump¬ 
tion currently has no alternatives: although conceptually 
problematic, it is seen as an operational necessity [i[23]- 
[25] . Short of resorting to individual-based modelling [4] , 
it remains unclear how the dynamics of ecological com¬ 
munities could be described in a naturally hierarchical 
way [iiaiiiiiiiiis]). 
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The hot debate surrounding the definition of the “unit” 
of ecological diversity prompts an intriguing question: is 
the partitioned community assumption an adequate de¬ 
scription of the microbial world? The species-based intu¬ 
ition is incontestably useful, but might it be forcing onto 
our data a structure that it does not possess [30||3T]? Ul¬ 
timately, of course, this question should be settled by ex¬ 
periments, and reports are conflicting [32H36] . However, 
ideally we should be asking not whether a species-based 
picture is “adequate”, but whether it is superior to alter¬ 
natives. Can we imagine an ecology where no grouping of 
individuals is privileged, and how would we describe it? 
Regardless of our stance on the “species problem”, we 
must recognize this question as an essential theoretical 
exercise: it is impossible to test a hypothesis unless we 
can imagine it to be false. Without an alternative, the 
partitioning assumption will continue to shape our ques¬ 
tions, our models and even our data analysis in the form 
of “operational taxonomic units”, an ill-defined proxy for 
a “bacterial species” extensively critiqued elsewhere m- 

Borrowing an idea from condensed matter physics, this 
work proposes a new way to tackle the species problem, 
approaching it “from the other side”, as illustrated in 
Fig.[T] Consider a heterogeneous community where some 
individuals are more alike than others, but the clustering 



FIG. 1. Clustering organisms into discrete species is often 
problematic (B). Rather than describing the natural world 
as a deviation from a perfectly clustered case (A), this work 
proposes a theoretical construct where structure can be grad¬ 
ually added to a fully disordered “ecology without species” 
(C). 
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FIG. 2. Coarse-graining degrees of freedom through 
clustering. A. The most common approach to coarse- 
graining the description of an ecosystem is by merging 
branches of a phylogenetic tree. B. An elementary cluster¬ 
ing operation where the abundances of two rabbit varieties 
are combined into a single degree of freedom. 


is imperfect, as is arguably the case in most natural situa¬ 
tions mnaEniEzi. Currently, the only existing approach 
treats this scenario as a small deviation from the clas¬ 
sic picture of well-clustered species. In physical terms, 
this is a perturbative expansion, where we start from 
a perfectly partitioned world (Fig. and attempt to 
recognize that species boundaries are “fuzzy” (Fig. m 
This work proposes an explicit construct where the per¬ 
turbative expansion can be performed around a different 
origin, starting from an unstructured phenotype back¬ 
ground (Fig. EP)- and adding some structure. 

Similar ideas have a long history and have been enor¬ 
mously productive in many areas of physics (cf. the clas¬ 
sic theory of metals that treats electrons as a free gas). 
The contribution of this work is to port this idea into an 
eco-evolutionary context. In order to achieve this, it is 
necessary to first develop a formalism that avoids postu¬ 
lating a community partitioning. This work extends the 
notion of an “ecomode” [38j |39] to construct a naturally 
hierarchical, rank-free description of community dynam¬ 
ics in terms of population fluctuation eigenmodes. This 
approach reduces to the familiar species-based perspec¬ 
tive when it applies, but remains well-defined even when 
the species description breaks down. The basic idea be¬ 
hind this formalism is that species clustering can be rein¬ 
terpreted as a special form of a more general operation, 
namely a change of basis in the compositional space. 


dom, such as abundance of individual strains, are as¬ 
sumed to be located at the leaves of a tree (based on phe¬ 
notypic similarity, phylogenetics, etc.). A coarse-grained 
description is then constructed by choosing a transect 
cutting across this tree, and combining branches lying 
downstream of the transect (Fig. |2}\). For concreteness, 
consider an example of combining just two leaves. Imag¬ 
ine a community where a species of rabbits comes in two 
varieties that differ only slightly: say, a “white” rabbit 
W and a “brown” rabbit B (Fig. [^). When describ¬ 
ing their predator-prey dynamics with foxes, a two- or 
three-variable description could be used, depending on 
the desired level of detail. The comprehensive three- 
variable description might capture come subtle features, 
but for many questions, the classic Lotka-Volterra model, 
whereby the two varieties are clustered together into a 
single variable “rabbits” i?, would provide a sensible ap¬ 
proximation: 


{W,B}^R. ( 1 ) 

This work proposes a different perspective on the clus¬ 
tering operation §■ To a first approximation, most en¬ 
vironmental perturbations will affect the population of 
rabbits as a whole; e.g. an increase in predator efficiency, 
or habitat deterioration will depress the population of 
both varieties together. Upon careful examination, how¬ 
ever, the small functional differences will manifest them¬ 
selves on a slow timescale: for example, if the predator is 
better attuned to spotting a particular coat colour, one 
variety may eventually outcompete the other. In this 
way, compositional changes naturally align not with the 
Lotka-Volterra axes W and 5, but with the rotated axes 
that provide an hierarchically-informed description: one 
variable R = W B (“rabbits”) that changes on a fast 
timescale, and another C = W — B (“coat color”) that 
changes on a slow timescale. This work proposes to rein¬ 
terpret Fig.[^ not as clustering, but as a change of basis, 
followed by the decision to ignore a slow-timescale vari¬ 
able. As we will see, an appropriately rotated basis can 
be constructed for an arbitrary population of interact¬ 
ing individuals, and can serve as a naturally hierarchical 
characterization of a community. Borrowing the term 
introduced in [38], one could call this the “ecomode de¬ 
scription” . 


I. REINTERPRETING SPECIES CLUSTERING 
AS A CHANGE OF BASIS 

Describing an ecosystem is a problem of dimensionality 
reduction. Given a community of billions of individuals, 
or millions of strains, the challenge is to find a reduced 
set of degrees of freedom that could allow characterizing 
its properties or dynamics with an acceptable degree of 
accuracy. 

The standard approach to this problem is through clus¬ 
tering. Most commonly, the microscopic degrees of free- 


II. THE EIGENMODE DESCRIPTION: 

THE GENERAL DEFINITION 

Let C denote a full individual-based description of a 
community: a list of all present organisms (labelled by 
indices /i, u...) and their functional characteristics. A 
general model for ecological dynamics can be written as a 
rule that associates to each individual /i a “rate of abun¬ 
dance change” r^: the probability per unit time with 
which it will either generate another individual (r^ > 0) 
or die (r^ < 0). In a most general model, this rate can 
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depend on any detail of the community state, as well as 
on external environmental parameters denoted by £\ 

ri, = r^{£,C) 

Within any such model, we can define a matrix Ai of 
interactions between individuals: is the effect that 

a hypothetical duplication of individual /i would have on 
the instantaneous division/death rate of u. 

= ru{£, {C with {i duplicated}) — ry{£,C) (2) 

Importantly, here fi and v are individuals^ not species, 
and so the definition of M. does not require the parti¬ 
tioned community assumption. For the classic Lotka- 
Volterra scenario, the “functional characteristics” of an 
individual need only specify whether it is a fox or a rab¬ 
bit, but we can allow this microscopic description to be 
arbitrarily detailed, down to the point of all individu¬ 
als being unique. Note that the interaction matrix A4 is 
defined instantaneously; the number of individuals and 
thus the size of A4 can change in time. 

The primary focus of this work is the eigensystem of 
the interaction matrix A4. For a population of J\f individ¬ 
uals, this Af X Af matrix will always have Af eigenvalues. 
Note, however, that if individuals fi and /i' happen to be 
functionally identical, then AA^y = Ai^'y for all and 
Ad has a zero eigenvalue. For a coarse model ascribing 
all individuals to only K ^ Af distinct phenotypes, the 
matrix Ad will be highly degenerate with only K non¬ 
trivial eigenvalues; the remaining Af — K will be strictly 
zero. A detailed functional description recognizing all 
individuals as unique will remove the degeneracy; how¬ 
ever, if the partitioning into K categories were indeed 
a good approximation, Af — K eigenvalues will remain 
small (here and everywhere, characterizing an eigenvalue 
as “small” refers to its absolute value). We are led to 
hypothesize that the spectrum of Ad may provide a nat¬ 
urally hierarchical description of ecological relationships 
in a community, with progressively smaller eigenvalues 
resolving finer and finer details. The grouping of individ¬ 
uals into discrete types, when it exists, is established as 
a result of this analysis; in this framework, it no longer 
constitutes a fundamental assumption. 


III. RESTRICTING GENERALITY: 
EIGENMODES IN LOTKA-VOLTERRA MODELS 

In the interest of full generality, the definition of the 
previous section was explicitly individual-based. The ex¬ 
istence of such a definition is a very important feature of 
the eigenmode framework; however, the inconvenience of 
dealing with discrete individuals is a technical nuisance 
that could obscure the larger points in the subsequent 
discussion. To build intuition about the eigenmode per¬ 
spective, the best approach is to illustrate it on simple 
examples. 


The most familiar setting is that of a generalized 
Lotka-Volterra model, with K continuous degrees of free¬ 
dom. Its linearized dynamics in the vicinity of an equi¬ 
librium point take the form: 

^ = J2MabSnk. (3) 

6=1 

Here Sn denotes the deviation from equilibrium abun¬ 
dance, and the index a runs from 1 to K. For simplicity, 
this is the context that will be used for most of the dis¬ 
cussion. Positing a finite number of K interacting phe¬ 
notypes may seem to contradict the stated goal of avoid¬ 
ing the “species” partitioning. However, the reader is 
invited to think of these phenotypes not as a “species” 
label, but as a microscopic, sub-species description; e.g. 
interacting bacterial strains that may or may not group 
into species. In this context, the variable-size individual 
interaction matrix Ad has only K non-trivial eigenvalues 
and eigenvectors, which correspond to the eigensystem 
of the Lotka-Volterra interaction matrix M (see the Sup¬ 
plementary Material; SM). The eigenvalues have units of 
inverse time and correspond to the relaxation time scales 
of perturbations applied along each eigenmode. 

In order to make the eigenmodes interpret able, rather 
than positing the interaction matrix M in some arbitrary 
manner, it is extremely helpful to construct it from a 
functional description of what these phenotypes do (re¬ 
call the example above, where the W ^ B axis could 
be conveniently interpreted as “rabbits” and IF — H as 
“coat color”). One simple way to do this is provided by 
the MacArthur model of resource competition [4Ql[4^ . 
where phenotypes are characterized by the resources they 
are able to consume. This defines a particular instance 
of a competitive Lotka-Volterra model, while also giving 
us an intuition for phenotypic similarity, and a functional 
interpretation for the eigenmodes we will observe. 

Specifically, following ref. [43], consider a well-mixed 
microbial community in a medium where a single limit¬ 
ing element (e.g. carbon) exists in N forms (“resources” 
i G {1.. .V}; in what follows N = 20). Each resource 
can be metabolized with a dedicated “pathway”. The 
phenotypes are defined by two pieces of information: 
first, their requirement for the limiting element, and sec¬ 
ond, the pathways they carry. The former will be ran¬ 
domly drawn for simplicity, but the latter will be chosen 
strategically, by hand, in order to construct the most in¬ 
structive examples. The set of K phenotypes defines a 
K-hy-N binary matrix a of pathway presence/absence 
{aai = 1 if phenotype a carries the pathway enabling 
it to consume resource i), and this matrix a determines 
the phenotype-phenotype interactions. Specifically, the 
competition between a pair of phenotypes (a, b) is the 
stronger, the larger the overlap in their resource prefer¬ 
ence '^^craiO'bi‘ Under MacArthur’s model, such com¬ 
petition always results in a unique and stable equili- 
birum mills]. The mathematical details are provided 
in the SM, but for the discussion that follows, this de¬ 
scription should already be sufficient. 
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FIG. 3. For a well-clustered set of phenotypes, the eigenmode representation captures the hierarchy of phenotype similarity. 
A: Pathway presence matrix for 12 phenotypes, generated as variations on 3 “core phenotypes” (colored lines); phenotypes 
carrying a given pathway (shown in yellow) are able to metabolize the respective resource. The tree represents hierarchical 
clustering based on Hamming distance. Dotted transect corresponds to the clear grouping into 3 clusters. Under MacArthur’s 
model, resource competition in this ecosystem leads to a unique stable equilibrium. B: The spectral structure of the ecosystem 
shown in A, in the vicinity of its equilibrium point. Left: the eigenvalues of the Jacobian. All eigenvalues are negative, 
corresponding to a stable equilibrium; shown is the absolute value. Right: the eigenvectors, ordered by decreasing |A|; each is 
a linear combination of 12 phenotypes and is difficult to interpret directly. C: The same eigenvectors, projected into pathway 
expression space. It is now clear that the three dominant modes span the three “core phenotypes” (colored lines as in panel 
A). The fourth mode (arrowhead) corresponds to a subdivision of the first cluster (compare with panel A). 


IV. THE EFFECTIVE NUMBER OF SPECIES IS 
DEFINED BY A SPECTRAL GAP 

For the first example, consider a community composed 
of phenotypes depicted in Fig. [^, competing for 20 re¬ 
sources supplied in equal abundance. These 12 pheno¬ 
types were purposefully generated to form three clear 
clusters (see SM for the exact procedure; Matlab scripts 
(Mathworks, Inc.) reproducing all examples and figures 
are available as Supplementary file 1). This community 
can be interpreted as harbouring three species; within 
each species, phenotypes share a core set of pathways 
(indicated) and differ only in a small subset, similar to 
the “core” and “accessory genome” of a species [44] . How 
is this scenario seen from the eigenmode perspective? 

Resource competition between these phenotypes de¬ 
fines a 12-by-12 interaction matrix. Its eigenmode spec¬ 
trum is presented in Fig. [^. Each eigenvector is a lin¬ 
ear combination of phenotypes, and is difficult to inter¬ 
pret. However, since each phenotype is defined by its 
resource consumption, their linear combination can be 
conveniently characterized in the same way: the pathway 
presence/absence matrix a can be used to project each 
eigenvector into the pathway expression space (Fig. IT)- 
This provides a functional interpretation for the eigen- 
modes: each row in the matrix shown in Fig.|^ indicates 
how the community-level pathway expression changes if 
its composition is perturbed in the direction of a given 
eigenmode. 

We observe that the spectrum is dominated by 3 
modes, whose projections in the pathway expression 
space are linear combinations of the “core genomes” of 
the three phenotype clusters (Fig. IT, top 3 rows). This 
indicates that at the short timescale, the response of the 
system to a perturbation can be well approximated us¬ 
ing only 3 degrees of freedom. Above, characterizing this 
habitat as “harbouring 3 species” was an interpretation 


that relied on a subjective judgement of Fig. IT by a 
human observer. The eigenmode spectrum contains this 
information, and more, in an objective format, decom¬ 
posing community dynamics over a range of time scales. 

The first eigenvector is a collective mode, with each 
entry close to 1. This corresponds to the statement that, 
if we insisted on reducing the community to a single de¬ 
gree of freedom, our best choice would be to say that all 
resources are being consumed at a certain overall rate. 
However, this is a poor approximation; capturing the 
dominant fast dynamics requires 3 degrees of freedom. 
Far below the dominating eigenvalues, a fourth one be¬ 
comes apparent (arrowhead). Examining the fourth row 
in Fig. and comparing it to panel A reveals that this 
mode corresponds to a finer structure within the first 
cluster, where coexisting variants differ in their ability to 
consume resources 1 and 2. The eigenmode decomposi¬ 
tion provides a way to describe dynamics with increasing 
level of detail; smaller differences manifest themselves as 
smaller eigenvalues. Thus, the eigenmode description of 
community structure is naturally hierarchical. 

The key property of the eigenmode spectrum of Fig.[^ 
is the existence of a noticeable gap in its spectrum, which 
provides a justifiable cutoff point for the number of de¬ 
grees of freedom to include into a low-dimensional de¬ 
scription (note that in some contexts, the term “spectral 
gap” is used in a narrower sense to describe the difference 
between specifically the first and second eigenvalues; this 
is not the meaning adopted here). Whenever an ecosys¬ 
tem exhibits a spectral gap, one could define the effective 
number of species as the number of dominant eigenval¬ 
ues. Importantly, however, the spectral gap need not 
be unique, and is a fluid property: under a continuous 
change of parameters, gaps can move, appear and disap¬ 
pear. The next example explores the implications of this 
fact. 
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FIG. 4. A smooth interpolation between two unrelated clustering patterns into, respectively, two and three effective species. A: 
Pathway presence matrix (same format as in Fig.[^) for 11 phenotypes. Their preference for resources 1-3 induces a grouping 
into three clusters, but the remaining pathways induce a conflicting clustering pattern (with only two groups). B: If the supply 
of resource 1-3 is low, community dynamics exhibits two dominant modes, reflecting the 2-cluster structure. As the supply 
of these resources is increased (panel C), the preference for resources 1-3 becomes the more relevant clustering criterion for 
describing the dynamics. D: The same 11 phenotypes now behave as a community of 3 effective species. Unlike a speciation 
event, where the uncertainty in defining the number of species is merely a matter of the desired level of detail (c/. Fig. [^, the 
intermediate regime of panel C interpolates between unrelated clustering patterns (not nested). 


V. MOVING THE SPECTRAL GAP: 

“TWO AND A HALF SPECIES”? 

Consider now a community of 11 phenotypes of 
Fig. |1\. For equiabundant resources, we observe the 
expected “2-species” spectrum (Fig. |^), the dominant 
eigenmodes reflecting the two “core genomes” that cor¬ 
respond to the dominating binary clustering pattern. 

However, by construction (see SM for details), each of 
these phenotypes specializes in only one of the first three 
resources. If the supply of these resources is strongly 
increased, this preference will become the more relevant 
clustering criterion; the same 11 phenotypes now behave 
as a community of 3 effective species (Fig.|^). Note that 
in this regime, the eigenmode corresponding to the binary 
clustering (which used to be mode #2) is again part of 
the spectrum (now as mode #4), but is characterized by a 
much smaller eigenvalue, signaling much slower dynamics 
than those governed by the dominant three resources. 

Both these extremes were easy to interpret, but how 
should we think of the regime lying in between (Fig. Ep)? 
A smooth process leading to a change in the number of 
species may seem familiar from the study of speciation. 
One cluster of phenotypes can be smoothly separated into 
two, and if we insisted on counting species, then some¬ 
where in between there must be a gray area of “more 
than one, but less than two”. This inability to draw a 
hard threshold of exactly when two close varieties should 
be considered distinct species is not, however, a signifi¬ 
cant challenge to the species paradigm. Whether a finer 
subdivision is warranted is merely a question of the de¬ 
sired level of detail. This was already illustrated in the 
previous example: the spectrum of Fig. exhibits two 
significant gaps, leaving us the choice of characterizing it 
with either three or four degrees of freedom. And just 
as in Fig. subjecting the community of Fig. Ik to an 
increased abundance of resources 1 and 2 would reduce 
the first spectral gap, making the finer (4-species) de¬ 


scription increasingly more justifiable. In this scenario, 
the origin of our inability to define an exact number of 
species is well understood. 

However, Fig. shows that the problem is more gen¬ 
eral. The somewhat artificial example of Fig. IIP was 
purposefully constructed to show that a smooth change 
of parameters can interpolate between two clustering pat¬ 
terns even when they are entirely unrelated (not nested): 
the same 11 phenotypes are grouped as {1-5}, {6-11} in 
one case (panel B), and {1, 6}, {2, 7-10}, {3-5, 11} in 
the other (panel D). The “gray area” in between is of a 
different nature than the uncertainty about whether to 
further subdivide a particular cluster. This suggests that 
the intermediate regime of Fig. HP should be interpreted 
not as “two and a half species”, but rather as a situa¬ 
tion where a gap-less spectrum means that the effective 
number of species cannot be defined at all. 

Characterizing such a “regime without species” is in¬ 
deed our ultimate destination. Unfortunately, this par¬ 
ticular example suffers from an obvious problem: it was 
constructed within a Lotka-Volterra framework with 11 
degrees of freedom. Whether or not they can be grouped 
into clusters, such a model can always be seen as describ¬ 
ing an interacting system of 11 species. 

This issue is a consequence of the decision, made early 
on, to work within the familiar Lotka-Volterra context. 
No simple Lotka-Volterra model can generate a truly gap¬ 
less spectrum: as discussed above, a model grouping in¬ 
dividuals into K discrete categories always has a gap sep¬ 
arating the K non-trivial eigenvalues from the vast num¬ 
ber of those that are strictly zero, as a consequence of 
the partitioning assumption. Similarly, virtually all clas¬ 
sical ecological models describe systems whose eigenvalue 
spectrum exhibits a prominent gap. What reasons do we 
have to believe that natural ecosystems, e.g. the micro¬ 
bial communities inhabiting our intestines, possess this 
property? Absent such evidence, we must acknowledge 
that spectral gaps, whose existence is a prerequisite for 
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FIG. 5. Constructing a model case of an ecology without species. A: A multi-patch extension of MacArthur’s resource 
competition model. Externally supplied resources are split between P patches harboring separate communities (here, P = 
3). Patch boundaries are assumed to be (weakly) permeable to resources, but not organisms (see SM for details). In this 
cartoon, phenotypes are represented by circles; size signifies abundance, while colors correspond to resources consumed. B: 
The phenotype-phenotype interaction matrix for a particular realization of the model with A = 20 resources and P = 15 
patches; the 298 coexisting phenotypes are ordered by patch. Organisms in the same patch interact more strongly than across 
patches, leading to a block-diagonal structure. C: The top P eigenmodes reflect the P-patch structure of the habitat. Modes 
2 through P have vanishingly small projection into the pathway expression space: they represent inter-patch dynamics and 
are hidden from an external observer who cannot resolve individual patches. D: The shape of the eigenvalue distribution 
is independent of P (left); modes “hidden” from external observer (2 through P) are omitted. The unique dominant mode 
is the collective mode (arrowhead). The remainder of the spectrum follows a dense distribution with no gap, reflecting the 
unstructured background of interacting phenotypes. 


a low-dimensional “species-based” description to be ade¬ 
quate, could be an artifact of familiar theoretical models, 
and not necessarily the best approximation of ecological 
reality. If one could build an ecosystem model whose 
spectrum really did exhibit no gaps, such a scenario fun¬ 
damentally could not be conceptualized as a system of in¬ 
teracting species. One signature of this regime is that any 
attempt to approximate it with a iF-dimensional Lotka- 
Volterra model will appear to contain exactly K species, 
for any K. The next section will construct one concrete 
example exhibiting this behavior. 

VI. AN ECOLOGY WITHOUT SPECIES 

Constructing an ecosystem with a dense spectrum re¬ 
quires a model where the number of distinct phenotypes 
can be arbitrarily large. It would be convenient to also 
maintain two useful features of the previous examples: 
a functional characterization of phenotypes, and the ex¬ 
istence of a stable equilibrium. Finally, to underscore 
the interpretation of this scenario as a regime “without 
species”, we will set an additional objective, requiring the 
set of coexisting phenotypes to be unstructured, so that 
no grouping is privileged. Ideally, then, we would like 
to have a model where an arbitrarily large, unstructured 
set of functionally defined phenotypes can co-exist at a 
stable equilibrium. 

In the classic (well-mixed) framework of MacArthur, N 
resources can stably sustain at most N phenotypes [42] . 
However, imagine a spatially structured habitat with P 
patches, whose boundaries are permeable to resources, 
but not to organisms (Fig.|^). The mathematical details 
of this model are provided in the supplement; however, 
qualitatively, the behavior is intuitive. At infinite perme¬ 
ability (resources are exchanged freely), the P patches are 


effectively a single patch, and we are back to the classic 
MacArthur model. If permeability is zero, the P isolated 
patches can sustain up to PN phenotypes, but these have 
no interactions across patch boundaries. At intermediate 
permeability, we find the desired situation where the in¬ 
teractions are nontrivial, and yet the total number of 
coexisting phenotypes can be arbitrarily large. Although 
the number of distinct binary vectors is bounded by 2^, 
two phenotypes located in different patches count as dis¬ 
tinct even if their pathway content is identical, because 
they interact differently with other phenotypes (for ex¬ 
ample, of the two “red-blue” phenotypes in the cartoon in 
Fig.|5}\, one has a strong competitor in the same patch, 
but the other does not). 

For each patch, define an independent pool of local 
phenotypes that is fully random (each phenotype carries 
or not a given pathway with probability 1/2). After the 
system equilibrates, the total number of coexisting phe¬ 
notypes is of order PN. Fig. shows the complete 
interaction matrix for 298 phenotypes in P = 15 patches 
(see SM for details). Predictably, phenotypes group by 
patch: at weak resource permeability, organisms interact 
more strongly within a patch than they do across patches. 
However, the pathway content of these phenotypes is, by 
construction, completely unstructured. 

These observations are reflected in the eigenmode 
structure. P eigenvalues clearly stand out in the spec¬ 
trum (Fig. [^), reflecting the patch structure: at the 
very top, we recognize the collective mode; the remaining 
P — 1 modes are pairwise differences between the collec¬ 
tive modes of the individual patches. These dominant 
modes are followed by a dense spectrum of eigenvalues, 
whose distribution exhibits no gap (Fig. |^): aside from 
the patch structure, no grouping of phenotypes is in any 
way privileged, offering no possibility for a coarse-grained 
description. 
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In order to take the limit P ^ oo, consider now the 
projection of the eigenvectors into the A/'-dimensional ex¬ 
pression space. This corresponds to adopting the per¬ 
spective of an observer who can measure the community¬ 
wide expression of metabolic pathways, but cannot re¬ 
solve individual patches. Such an observer can see 
the dominant collective mode, but modes 2 through P 
are effectively hidden: their projection into the path¬ 
way expression space has vanishingly small components 
(panel C). To understand this, consider a perturbation 
where all phenotypes in patch 1 experience an increase 
in abundance, and all phenotypes in patch 2, a decrease. 
Such a perturbation will relax quickly, corresponding 
to one of the dominant (strongly negative) eigenvalues. 
However, the internal dynamics represented by this mode 
is effectively hidden from an outside observer whose mea¬ 
surements do not resolve individual patches. Such an ob¬ 
server can only see a unique collective mode over a fully 
unstructured background (modes from P 1 onwards): 
a dense spectrum, whose shape is independent of P. 

This characterization allows sending P to infinity. This 
limit is proposed as a tractable model case of an “ecology 
without species”. In this regime, the number of interact¬ 
ing phenotypes becomes infinite, their pathway content 
is unstructured, and the spectrum is dense; yet this con¬ 
struct continues to describe a valid ecology at equilibrium 
whose properties could be investigated. For instance, if 
the supply of some resource is increased, this system will 
have some well-defined response, which our hypothetical 
observer could measure and could attempt to predict. 

What is the utility of this regime? The extreme, fully 
unstructured case is certainly a poor approximation of 
any natural setting. However, it provides a baseline that 
can be perturbed to describe novel scenarios of “weakly 
structured” ecologies, as illustrated in the next and final 
example. 


VII. DESCRIBING WEAKLY STRUCTURED 
ECOLOGIES 

Fig. [^repeats the analysis of Fig. with one exception: 
the first and last pathway are forbidden from ever ap¬ 
pearing within the same organism. Specifically, for each 
patch, the procedure generating the local phenotype pool 
was modified as follows: each entry in the pathway pres¬ 
ence matrix a was independently set to 1 with probability 
1/2 as before, but whenever this resulted in a phenotype 
carrying the forbidden pathway pair, this phenotype was 
removed from the pool. This kind of structure is not im¬ 
plausible: a detrimental cross-talk between certain path¬ 
ways can make their simultaneous activity in the same 
organism difficult; for example, the enzyme nitrogenase 
is inactivated by oxygen and thus must be kept separate 
from oxygen-utilizing pathways. 

The spectrum of this ecosystem, projected into the 
pathway expression space, is shown in Fig. 1^. Com¬ 
paring to Fig. [^, we see that there are now two modes 


A 

1 

0.1 

0.01 



Pathways 


- P = 20 


- P = 50 

■ ■ 


-P= 100 


Eigenmode #P+1 


FIG. 6. A weakly structured ecology as a perturbation of the 
no-species regime. A: When the phenotype pool is modified 
to introduce an anticorrelation between first and last path¬ 
way, the spectrum acquires a mode reflecting this structure 
(highlighted). This new mode is the dominant observable 
mode (after the collective mode). C: The leading observable 
eigenmodes remain at the same position in the spectrum (ar- 
rrowheads) and retain their structure (right). The rest of the 
spectrum follows the same distribution as in the unstructured 
regime (shaded; compare with Fig.[^). 


that stand out from the gap-less background. In addi¬ 
tion to mode #1 (the collective mode), there is now also 
mode #(P + 1) that represents the only perturbation 
imposed on an otherwise unstructured phenotype back¬ 
ground. The spectrum in Fig. was computed for one 
particular random realization of a system with P = 15 
patches; however, as P is increased, the projection of 
this mode into pathway space and its position in the 
spectrum both remain invariant (Fig. |^). This exam¬ 
ple demonstrates that the fully unstructured regime con¬ 
structed here can indeed be meaningfully perturbed to 
describe the dynamics of “weakly structured” ecologies, 
an interesting theoretical scenario that was previously 
out of reach. 

To maintain consistency with the previous exam¬ 
ples, the approach taken here focused exclusively on re¬ 
source competition. The “hidden” modes correspond¬ 
ing to unobservable inter-patch dynamics is the price to 
pay for this simplicity. Other approaches are possible; 
one might, for instance, imagine constructing a dense- 
spectrum model by means of a fully individual-based de¬ 
scription, or through a sequence of iterative subdivisions 
of a Lotka-Volterra model (similar to the approach taken 
in Ref. [46], but repeated an infinite number of times). 
One of these methods might yield a more elegant con¬ 
struction; however, the multi-patch example constructed 
here has the advantage of being easily interpretable, and 
potentially even analytically tractable with methods de¬ 
veloped in Ref. [47] . 


VIII. DISCUSSION 

The mathematics of the “eigenmode approach” de¬ 
scribed here is very simple. Anyone familiar with the 
Principle Component Analysis (PCA) has encountered 
the problem of choosing the number of components on 
which to focus. In rare circumstances, a few components 
















are sufficient to capture most of the variance; the generic 
situation, however, is that the spectrum shows no obvi¬ 
ous gap, offering no natural truncation point [48]. Nor¬ 
mally, one thinks of the PC A modes as linear combina¬ 
tions of the fundamental degrees of freedom, namely the 
species abundances. Here, the novelty is to ask whether 
the eigenmodes could be used as the fundamental degrees 
of freedom themselves. 

Even if the concept of species is not called into ques¬ 
tion, the eigenmodes provide a naturally hierarchical de¬ 
scription which allows asking new questions. Are there 
general laws governing the structure and evolution of 
community spectra? What ecological mechanisms might 
cause a given eigenvector to change in a predictable way, 
or a cause a given eigenvalue to become more prominent? 
How is the spectrum affected by evolutionary mecha¬ 
nisms? One may, for example, expect that purifying 
selection encourages the development of spectral gaps, 
while diversity-creating mutations have the opposite ef¬ 
fect. A major advance in classical physics came with the 
development of the renormalization group, a framework 
explaining why many physical systems, as they become 
large, do in fact develop a gap in the spectrum, enabling 
some low-dimensional “effective theory” to become an 
excellent description [49] . Ecology and evolution provide 
a new context for asking similar questions, which are par¬ 
ticularly relevant today with the study of highly diverse 
microbial ecosystems of medical and environmental in¬ 
terest nsiiTTj. 

The eigenmode formalism admits the possibility that, 
for a given ecosystem, the most relevant degrees of free¬ 
dom may take a different form than simply the combined 
abundance of some taxa. As a simplest example, imag¬ 
ine a two-species consortium Si, S 2 that runs the same 
metabolic reaction that species 83 performs alone. In 
this scenario, the combination Si ^ S 2 — S 3 might be a 
functionally “negligible” (slow) eigenmode, but the re¬ 
duction of dimensionality this observation affords does 
note reduce to simple branch-merging on a phylogenetic 
tree. Situations when one diverse set of species is essen¬ 
tially interchangeable with another diverse set are likely 
the rule, rather than an exception, as evidenced by the 
dramatic compositional variations of human-associated 
microbiota across healthy individuals, at all taxonomic 
levels [16]. This suggests that the number of negligible 
eigenmodes may indeed be large, and recognizing them 
could simplify the task of modelling these communities 
even within the traditional species-based framework. 

However, more profoundly, the perspective developed 
here allows conceptualizing a novel scenario of an “ecol¬ 
ogy without species”. The classic Lotka-Volterra model 
investigates an ecology of perfectly defined species, a 
regime that has been studied extensively. This work pre¬ 
sented a theoretical construct that can be seen as the op¬ 
posite extreme, one where the phenotypes are completely 


unstructured, and yet form a valid ecology. 

In any real setting, organisms sharing a recent common 
ancestor are functionally similar, and the approximation 
“every single organism is unlike any other” is even more 
extreme than the converse assumption of perfect clus¬ 
tering into a small number of species. Nevertheless, a 
theoretical investigation of this unstructured regime will 
likely be a valuable step towards understanding the mid¬ 
dle ground occupied by the real-life communities; in phys¬ 
ical terms, its role would be to provide a different origin 
for a perturbative expansion (Eig. [^. Eor example, in 
condensed matter physics it was recently proposed that 
a granular medium at the brink of a jamming transition 
can be thought of as a “maximally disordered” solid, the 
opposite extreme of a perfectly ordered crystal m- Im¬ 
portantly, it was shown that certain properties of real- 
life solids, even at relatively low disorder, can be more 
adequately described as a perturbation of this “maxi¬ 
mally disordered” regime, rather than as a perturbation 
of a perfect crystal m- Similarly, it is intriguing to hy¬ 
pothesize that some behaviors of highly diverse natural 
ecosystems may be better described as a perturbation of 
a fully unstructured no-species regime (adding some ten¬ 
dency of phenotypes to cluster), rather than as a small 
deviation from the classic picture of well-clustered species 
(where we attempt to recognize that species boundaries 
are “fuzzy”). This approach would provide a system¬ 
atic framework to investigate the implications of “struc¬ 
tured” and “unstructured” variation imposed at a range 
of scales, effects whose importance at the intraspecffic 
level was recently highlighted m- 

The eigenmode formalism is by no means limited to 
resource competition models; MacArthur’s model merely 
provided a convenient setting for constructing easily in¬ 
terpret able examples. The framework described here 
does, however, suffer from a number of other limita¬ 
tions: for example, only deterministic dynamics were 
considered, and the discussion was restricted to small 
fluctuations around some ecological state. Nevertheless, 
the question raised here was whether our description of 
ecosystems could move beyond the inherently discrete 
concept of a “species”. In light of this challenge, the 
simplifications above provide a reasonable starting point. 
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SUPPLEMENTARY MATERIAL 


Appendix A: The resource competition model 
(single patch) 

Examples in Figs. and were constructed within the 
resource competition framework used in Refs. [43l ? ]. 

For a detailed discussion of the metagenome partitioning 
model, its relation to the classic model of MacArthur im, 
and the proof of existence, uniqueness and stability of the 
equilibrium state of its dynamics, the reader is referred 
to Ref. [43] (main text and supplementary material). 
The procedure for efficiently computing this equilibrium 
state numerically, as a convex A-dimensional optimiza¬ 
tion problem, is described in Ref. [? ]. All computa¬ 
tions were performed in MatLab (Mathworks, Inc.), and 
scripts reproducing all examples and figures are available 
as Supplementary File 1. 

For the sake of completeness, the model will be briefly 
defined de novo. The subsequent sections will describe 
the specifics of how examples in Figs. and were 
generated. The extension to a multi-patch habitat and 
Figs. ii will be discussed in a separate section. 


1. Definition of the model 


Consider a microbial community in a habitat where 
a single limiting element X (e.g. carbon) exists in N 
forms (“resources” i G {1...A}) denoted A, 5, etc. 
The substrates can be utilized with “pathways” Pi (one 
specialized pathway per substrate). A phenotype a is 
defined by the pathways that it carries. Specifically, its 
“pathway content” is a binary vector of pathway pres¬ 
ence/absence: aa = {(Jai} = {1, 1, 0, 1, . . . }. 

For a given community, let Ti be the total number of 
individuals capable of utilizing substrate i: 






Here Ua is the abundance of phenotype a in the commu¬ 
nity (number of individuals with this phenotype). As¬ 
sume a well-mixed environment, so that each of these Ti 
individuals gets an equal share Ri/Ti of the total ben¬ 
efit Ri (e.g. carbon content) available from resource i. 
Any one resource is capable of sustaining growth, but 
accessing multiple cumulates the benefits. The fate of an 
individual of phenotype a is determined by the resource 
surplus Aa it experiences: 


Aa = - Xa- 


(SI) 


Here the first term is the benefit harvested by all carried 
pathways, and Xa represents the metabolic “cost” of phe¬ 
notype a, namely its requirement for the growth-limiting 
element A. The resource surplus A is “spent” on gener¬ 
ating new biomass. Equating, as in Ref. [43], the biomass 


of an organism with its cost for simplicity, we posit that 
for an individual of phenotype a, its division/death rate 
(as defined in the main text) is given by: 

ra = —Aa, (S2) 

^OXa 

where tq sets the unit of time. In terms of phenotype 
abundances, this dynamics translates into: 

^ = ffa(n) = naAa. (S3) 

at ToXa 

Here n is the vector of abundances of all phenotypes, 
and the function ga{n) is defined by the right-hand side 
of this equation. 

One does not need to make the particular choice of 
equating phenotype biomass with its metabolic cost; 
other choices are possible. The same applies to the choice 
of resource competition as the unique ecological interac¬ 
tion on which to focus. Here and elsewhere, no claim is 
made that the particular choice made here is the best 
representation of biological reality. As explained in the 
main text, the examples used in this work were generated 
specifically to illustrate certain features of the eigenmode 
perspective, and most notably the role played by the gap 
in the model’s spectrum. Many simplifying choices are 
made in order to make these illustrations easiest to in¬ 
terpret. 


2. Stability analysis 


The dynamics (S3) always lead to a unique, stable equi¬ 
librium m? ]. This equilibrium may consist of fewer 
phenotypes, as some phenotypes may go extinct during 
equilibration. Denote S the set of surviving phenotypes; 
for a G S' the equilibrium condition reads > 0, = 0. 

For the phenotypes that went extinct (a ^ S), we have 

na = 0, Aa < 0. 

We can now perform the stability analysis of the dy¬ 
namics (S3) around the equilibrium point, i.e. consider 


the eigenmodes of the matrix 


Mab = ^. (S4) 

dub 

Considering separately the types that are present and 
those that are absent, one can write: 


Mab = 


Ur 


Ri 


1 


^ ^ ^ai^hi rp2 if ^ 
a • -^i 


(S5) 


roXa 


'^ab^a 


\i a d S. 


Since the equilibrium is stable [43l ? ], all eigenvalues 

of this Jacobian matrix are negative. Restricted to the 
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space of phenotypes that are present, this matrix is given 
by: 

^abla^beS ~ ~ (S6) 

mXa . J-i 

Our object of interest is the eigensystem of this matrix. 

3. Lotka-Volterra vs. an individual-based 
description 

The main text considered two approaches to describ¬ 
ing ecosystem dynamics: the discussion began with a 
general individual-based description (with an individual- 
individual interaction matrix Al), which was then re¬ 
stricted to Lotka-Volterra-type dynamics for simplicity 
(with a phenotype-phenotype interaction matrix M). 
The main text claimed that if individuals fall into just 
K categories, then the non-trivial eigenvalues of the 
individual-individual interaction matrix M. in the vicin¬ 
ity of equilibrium coincide with the eigenvalues of the 
K-hy-K matrix M of the corresponding Lotka-Volterra 
model. (Assuming the abundances are sufficiently large 
that the continuous Lotka-Volterra model is a valid ap¬ 
proximation.) 

This claim is easy to verify. As above, let Ua denote 
the abundances of the K phenotypes, and let Mab be the 
matrix that defines linearized dynamics in the vicinity of 
this equilibrium point: 

^^ = y2MabSnb. (S7) 

6=1 

Here Sn is the deviation from equilibrium abundance 
n^a^. To convert this into an individual-based descrip¬ 
tion, we “split” each phenotype into ria individuals, con¬ 
sidered separately. What is the corresponding individual- 
individual interaction matrix? For two individuals, /i of 
phenotype a and u of phenotype 6, the interaction 
is, by definition, the change in growth/death rate of u 
that would be caused by a duplication of individual ji. 
Let us make two observations: first, duplicating fi means 
increasing of by 1 unit. Second, the population growth 
rate of a phenotype b and the division/death rate of an 
individual u are related by the factor rij^. Therefore: 

= -Mab. (S8) 

Ub 

The final observation is that each eigenvector Va of the 
matrix Mab can be converted into an eigenvector of the 
matrix by duplicating each entry ria times. The 

relation ( |S8| ) guarantees that the result will indeed be an 
eigenvector of M with the same eigenvalue. 

4. Figure 2 

The purpose of the example shown in Fig. 2 was to il¬ 
lustrate that when phenotypes exhibit a clear structure. 


the eigenmode formalism successfully captures both the 
broad clustering pattern, and the more subtle subdivi¬ 
sion. Accordingly, 4 “core phenotypes” were chosen, two 
of which are very close. The colors refer to the colored 
lines indicating the “core phenotypes” in Fig. 

1. carrying pathways 2 through 8 (green, variant 1); 

2. carrying pathways 1 through 8, but not 2 (green, 
variant 2); 

3. carrying pathways 9 through 14 (magenta). 

4. carrying pathways 15 through 20 (cyan); 

The set of 20 competitors was constructed by taking five 
copies of each core phenotype, and randomly flipping 
pathway presence/absence with probability 0.05 to gen¬ 
erate some diversity. Here, “flipping a pathway” means 
adding it to a phenotype where it was absent, or remov¬ 
ing if it was present. The cost of each competitor was 
then assigned as follows: 

Xa ~ ^ ^ ^ia ^Ca- 
i 

Here is a standard Gaussian random variable, and e 
was set to 0.01. In other words, all competitors have 
the same cost per carried pathway, with a small pertur¬ 
bation. A cost that scales linearly with the number of 
carried pathways ensures that phenotypes carrying more 
pathways do not automatically have a competitive ad¬ 
vantage [43]. 

The community of 20 competitors that resulted from 
following the above procedure (for some particular seed 
of the pseudo-random number generator) was then equili¬ 
brated (with all resources supplied at equal rate, = 1), 
following the numerical procedure described in Ref. [? ]. 
As a result of this equilibration, some phenotypes went 
extinct; the 12 survivors are shown in Fig. 

5. Figure 3 

The procedure generating this example is very simi¬ 
lar: first, a pool of competitors with desired properties 
is constructed by hand. The competitors are assigned 
costs with a slight random component to avoid any ac¬ 
cidental symmetries (e.g. two organisms whose costs are 
exactly the same). The community is equilibrated, and 
the pathway content of survivors is reported. 

For Figure 3, phenotypes were constructed as varia¬ 
tions on 6 “core phenotypes”. Their structure (put in by 
hand) can be read off from Fig. 3A: 

1. pathways 1 and 4-12; 

2. pathways 2 and 4-12; 

3. pathways 3 and 4-12; 

4. pathways 1 and 13-21; 
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5. pathways 2 and 13-21; 

6. pathways 3 and 13-21. 

Note that in this case N = 21 pathways were used; this 
was done to ensure that all six core phenotypes carried 
the same number of pathways. 

As before, a pool of 30 competitors was gener¬ 
ated by taking 5 copies of each core phenotype, and 
adding/removing a few pathways to generate some diver¬ 
sity. To preserve the carefully prepared structure, path¬ 
ways 1 through 3 were left intact, ensuring that all phe¬ 
notypes specialize in exactly one of these first three re¬ 
sources. Of the remaining 18 pathways (numbered 4-21), 
exactly two were randomly “flipped” for each phenotype. 

Finally, the pool of competitors was equilibrated for a 
range of resource conditions. The supply of resources 4 
through 21 was fixed at 1, while the supply of resources 1, 
2 and 3 was varied (synchronously) from 0.5 to 50. Only 
phenotypes that survived at non-zero abundance for the 
entire range of these conditions were retained. This pro¬ 
cedure ensures that throughout the range of conditions 
shown in Fig. 3C, the community consists of the exact 
same set of phenotypes, in this case 11 of them. 

At the end of this protocol, the resulting set of pheno¬ 
types was verified to contain at least one representative 
from each of the 6 core phenotypes. If this were not the 
case, the procedure would have been repeated for another 
choice of the seed for the pseudo-random generator. 

Appendix B: The multi-patch model: Figs and 
1. Definition of the model 

In the single-patch model above, one can define the 
availability of a resource as follows: 



Recalling the definition of the “resource surplus” 

^ ^ ^ai^i Xa: 
i 

we see that hi determines the benefit that any organism 
receives from carrying the relevant pathway. Consider, 
for example, a specialist that only carries pathway #1, 
and has cost xi- Introduced in a community where the 
availability of resource 1 is hi, this organism will be able 
to grow if and only if xi < hi. If this condition is sat¬ 
isfied, then this organism will multiply until it depletes 
resource 1 to exactly hi = Xi- We see that xi plays the 
role of Tilman’s R*, and hi can be interpreted as the 
concentration of resource i in the growth medium (up 
to some constant factor which ensures the right dimen¬ 
sion). This argument is described in detail in Ref. [? ]; 
the reason it is briefly reproduced here is to motivate the 
multi-patch construction that follows. 


Consider two communities in two isolated patches; one 
depletes resource i to the other to hp\ For con¬ 
creteness, let hp^ < h-^^ the second community is more 
efficient at depleting resource i, and drives its equilibrium 
concentration to a lower value. If the patch boundary 
were permeable to this resource, we would expect a net 
flux of resource from the patch with higher concentra¬ 
tion to a patch with lower concentration. In the simplest 
model, one can postulate that the magnitude of this flux 
(from the first patch into the second) is 

5Ri=p{hf^ 

where p measures the permeability of patch boundaries. 
Let 7? denote the global supply of resources (the vector 
notation is used to avoid explicitly specifying the index 
i, i.e. 7? = {hZi}). In isolation, both patches would 
receive an equal share of resource influx However, 
the permeability of patch boundaries causes resources to 
be redistributed: the first patch now receives — SR, 

and the second receives + SR. 

Generalizing this to P patches, we obtain the following 
model for resource consumption in a multi-patch habitat. 
For each patch, the ecological dynamics remain identical 
to those of a single-patch MacArthur’s model described 
above. The coupling between patches manifests itself 
only in the magnitudes of resource influx experienced 
by each patch. The global external supply 7? is evenly 
split between all P patches; however, each patch may 
experience additional influx or outflux depending on the 
resource uptake efficiency of its local community. 

Specifically, if a labels patches, define h* as the average 
resource availability across all patches: 



Patches where the availability (concentration) of resource 
i is lower than h* experience a net influx of this resource 
from other patches, where the availability is higher. The 
total resource influx experienced by a given patch is 
therefore: 

where 

5R^°‘'’ = p{hf^ -h*). 

At zero permeability (p = 0), the patch dynamics are 
fully independent. At infinite permeability (p ^ oo), any 
difference in local resource availability would generate an 
infinitely strong flux, eliminating such difference; in this 
limit, the availability of resources in all patches must 
be identical = h*, and the standard single-patch 
MacArthur model is recovered. Both these limits have 
the expected behavior. 
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Numerically, for weak permeability p, the equilibrium 
of this dynamics can be found through a simple iterative 
procedure: each patch is separately equilibrated at fixed 
resource supply, at which point the resource exchange 
fluxes are computed, the amount of resource available 
at each patch is updated, and the procedure is repeated 
until convergence. 

2. Figs and 

To generate the example in Fig. at each patch, a 
local pool of random 200 competitors for 20 resources 
was generated (at each patch, the pathway presence ma¬ 
trix was a random binary matrix of size 200-by-20, each 
entry set to 0 or 1 with probability 1/2). The number 
of patches was P = 15 for the example shown in panels 


B and C, and increased to 20, 50 and 100 to compute 
the eigenvalue distributions in panel D. As before, the 
cost of each phenotype was set equal to the number of 
pathways it carried, with a small random contribution 
to break any accidental symmetries (adding a Gaussian 
random variable with zero mean and width 10“^). The 
permeability parameter was set to 1 (a value empirically 
determined to have some weak effect on the set of co¬ 
existing phenotypes, inducing some interaction between 
patches), and the iterative equilibration procedure was 
performed until the relative magnitude of change in re¬ 
source influx experienced by each patch would fall below 
10 -^. 

To generate Fig. the exact same procedure was fol¬ 
lowed, with one exception: prior to the equilibration pro¬ 
cedure, any phenotype containing the “forbidden” path¬ 
way pair was removed from the competitor pool. 


